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Abstract 

In this paper we investigate the problem of a long self-avoiding polymer chain immersed in a 
random medium. We find that in the limit of a very long chain and when the self-avoiding inter- 
action is weak, the conformation of the chain consists of many "blobs" with connecting segments. 
The blobs are sections of the molecule curled up in regions of low potential in the case of a Gaus- 
sian distributed random potential or in regions of relatively low density of obstacles in the case of 
randomly distributed hard obstacles. We find that as the strength of the self- avoiding interaction 
is increased the chain undergoes a delocalization transition in the sense that the appropriate free 
energy per monomer is no longer negative. The chain is then no longer bound to a particular 
location in the medium but can easily wander around under the influence of a small perturbation. 
For a localized chain we estimate quantitatively the expected number of monomers in the "blobs" 
and in the connecting segments. 

PACS numbers: 36.20.Ey, 05.40-a, 75.10.Nr, 64.60.Cn 



I. INTRODUCTION 



The properties of polymer molecules trapped in a random environment like a gel network 
|l|] or inside porous materials or membranes 0, ^, ^ is of considerable interest. On the 
theoretical side many papers were written on the subject in an effort to understand the 
basic properties of the model, like the radius of gyration (or end-to-end distance) of a single 
molecule immersed in the random medium [5-14]. In several papers an ideal (Gaussian) 
chain has been used, which corresponds approximately to the experimental situation at the 
so called B-temperature when the solvent effectively screens the self-avoiding interaction of 
the chain. Even in that case the properties of a single molecule in a random potential are not 
simple. First, the properties depend on the type of randomness introduced- like a Gaussian 
random potential versus a sea of hard obstacles. Second, it may depend on whether the 
random potential is annealed or quenched, i.e. if the impurities are free to move or fixed. 
Additionally, different results apply if the chain is anchored or free to move. It turns out 
that for quenched randomness and a chain that is free to move, the size of the total volume 
of the medium is important for determining the chain size. For example, in the case of 
a random potential with a Gaussian distribution, the chain size saturates as a function of 
its length L (number of monomers) for a very long chain, in contrast to a dependence like 
Rg = aL^/'^ (where Rg is the radius of gyration and a is the monomer size) in the absence 
of a random potential . A logarithmic dependence of the chain size on the volume of the 
medium was argued in reference [0 using qualitative arguments. An analytic derivation of 
the logarithmic dependence was derived in reference using the replica method and the 



variational approximation, with a replica-symmetry-breaking solution. In reference |T3[ the 
result was justified from a mapping of the polymer problem to the problem of the localization 
of a quantum particle in a random potential. In reference results were obtained for a 
polymer immersed in a sea of random obstacles, and various behaviors were identified as a 
function of the system size. Again the use of localization theory and extreme value statistics 
were instrumental to the derivation. 

In some of the papers [|, |, || there was an attempt to include the effects of the self- 
avoiding interaction of the polymer. These attempts were far from complete. For example 
in Ref. it was assumed that the conformation of the polymer consists of one spherical 
blob and it was argued that a quenched random potential is irrelevant for a very long chain. 
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In Ref. analytical results were obtained for annealed disorder, and simulations were 
performed for strictly self-avoiding walks. Ref. presents numerical evidence for a size 
transition of the polymer as a function of the relative strength of the disorder and the self- 
avoiding interaction. The simulations were carried out for a random distribution of hard 
obstacles with a concentration exceeding the percolation threshold. It is the aim of the 
present paper to shed more light on this important problem. We will make use mainly of 
Flory-type arguments, and consider both the case of a Gaussian random potential and the 
case of randomly placed obstacles. A polymer with self avoiding interaction can not be 
mapped into a quantum particle at a finite temperature in a simple manner, because for 
a quantum particle there no impediment to return at a later time (or Trotter time) to a 
position it visited previously. 

An important point to keep in mind is the strength of the excluded volume interaction. 
If one consider a strictly self-avoiding walk on a lattice (SAW) corresponding to a non-self- 
intersecting chain, then the strength of the Edwards parameter v is fixed at 0(1) x 
where a is the step size (or monomer size) and depends only on the type of lattice. On the 
other hand one can consider a Domb- Joyce model [1^| where there is a finite penalty for 
self overlapping of polymer segments, and then the strength of v can be varied substantially 
and reduced continuously to zero. The interplay between the strength of the self avoiding 
interaction and the strength of the disorder can then be investigated to a larger extent. 
Experimentally the Edwards parameter is given approximately [rS| by f = a^(l — 2%), 
where x is the Flory interaction parameter, which depends on the chemical properties of the 
polymer and the solvent, and on the temperature (and pressure). It takes the value 1/2 at 
the G-point. The case x = corresponds to a solvent that is very similar to the monomer. 
In general good solvents have low x whereas poor solvents have high x resulting in v being 
negative. In the following we will restrict ourselves to the case of positive which leads to 
the more interesting and non-trivial results. 

Before we proceed to investigate the effects of the self-avoiding interaction, let us sum- 
marize in more detail the known results for a free Gaussian chain in a quenched random 
potential with a Gaussian distribution with variance and for a free Gaussian chain in a 
sea of fixed random obstacles with average concentration x per site. 

For an uncorrelated Gaussian random potential, it was argued using qualitative arguments 
in Refs. |]^ ^ that a very long Gaussian chain, that is free to move in the medium, will 
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typically curl up in some small region of low average potential. The polymer chain is said 
to be localized, and for long chains the end-to-end distance {Rf) becomes independent of 
chain length L (number of monomers) and scales like 

Rpoc{g\nV)-'^^'-''\ (1) 

where V is the volume of the system, d is the number of spatial dimensions, and the random 
potential satisfies (f/(x)f/(x')) = g6{'x — x'). The binding energy per monomer of the chain 
is given approximately by 

U^ind/Lr^iglnVf/^'-'l (2) 



These results were also obtained by the replica method in Ref. [12|, and rederived using a 



mapping to a quantum particle's localization in Ref. |T3[. For very short chains the end- 
to-end distance scales diffusively {Rp ~ L), and it saturates at the Rp value quoted above 
for large L. Notice, that in the infinite volume limit, the chain completely collapses. This 
results from the fact that the depth of the potential is unbounded from below, and the chain 
is always able to find with reasonable probability a deep enough narrow potential well to 
occupy, overcoming its tendency to swell due to the entropy of confinement. The collapse 
of the chain in the infinite volume limit agrees with the results for a chain in an annealed 
potential, since the ability of a free chain to scan all space for a favorable environment is 
equivalent to the random potential adapting itself to the chain configuration. 

In contrast, a tethered chain that is anchored at one end behaves very differently. Such 
a chain has a tadpole structure. The end of its tail is anchored and its head is a curled coil 
situated at a deep potential well. Such a chain has an end-to-end distance of order 

gl/(4-rf) I 

which is ballistic in three dimensions, since the chain has to look far to find a low potential. 
Here there is no dependence on the system's volume and the quenched and annealed cases 
are very different from each other even in the infinite volume limit. 



In reference we considered the effect of randomly distributed obstacles on the behavior 
of a free Gaussian chain. We have shown that in the presence of infinitely strong obstacles, 
that totally exclude the chain from visiting sites occupied by obstacles, with average concen- 
tration X per site, there are three possible behaviors of the end-to-end distance depending 



on the system's volume V. For 2 < c? < 4, if the volume of the system is smaller than 
Vi ^ exp(cia;~*^'^~^'*/^), then the polymer is localized, as in the case of a Gaussian random 
potential, and Ri ~ {x\nV)~^'^^~'^\ For Vi < V < V2, where V2 ~ exp (02x2/(^^+2) ) , 
the polymer size is given by Ru ~ (x/lnV)^^^'^ (ci and C2 are constants of order unity). 
Finally, for V > V2, the polymer behaves the same way as for an annealed potential, and 
RiH ~ {L/xf'^'^^^\ Resuhs are given to leading order in x for small x. These results 
are valid only when the average volume fraction of the obstacles {x) is smaller than the 
percolation threshold. When x is bigger than the percolation threshold we expect that the 
system breaks into independent domains whose volume is independent of, and generally 
much smaller than the total volume of the system. 

We end the introduction by revisiting the meaning of localization for polymers in a 
random medium. Although some authors connect the compact size of the chain when L ^ 00 
with the notion of localization, this is actually not so. The compact size should be viewed 
as a separate feature from the notion of localization. It is rather the binding energy of a 
chain Ubind that has to exceed the translational entropy InV. From Eq. (0) this amounts to 
the condition 

InV < L(^lnV)2/(4-'^), (4) 

which holds for any 2 < d < 4 when L is large enough (for 2 < d < 4 and any fixed L, 
the condition can be satisfied for large enough V) . This condition assures that the polymer 
will stay confined at a given location and will not, under a some small perturbation move 
to a different location. Thus repeating an experiment or a simulation with the same fixed 
realization of the disorder, but with different initial conditions, will result in finding the 
polymer situated at the same region of the sample as in a previous experiment, provided of 
course one waits enough time (which can be enormous) for the system to reach equilibrium. 
We observe that this condition is satisfied for large enough L provided the binding energy per 
monomer is positive. Another interpretation of the inequality given above in the context of 
equilibrium statistical mechanics is that the partition sum is dominated by the term involving 
the ground state as opposed to the contribution of the multitude of positive energy extended 
states. The contribution of these states is proportional to the volume of the system and thus 
the inequality above results from the condition 

exp(-L£'o) > V. (5) 
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What we will see in the following sections is that in the presence of a self-avoiding interaction, 
a localization-delocalization transition occurs when varying the strength of the the self- 
avoiding interaction for a fixed amount of disorder or alternatively upon varying the strength 
of the disorder for a fixed value of the self-avoiding interaction. 

II. CHAINS WITH SELF-AVOIDING INTERACTION-THE CASE OF A GAUS- 
SIAN RANDOM POTENTIAL 

We now study the case of a self-avoiding chain. In this section we consider the case 
of a random potential with a Gaussian distribution, and postpone to the next section the 
investigation of the case of a random distribution of hard obstacles. For simplicity, the 
discussion in the rest of the paper will be limited to three spatial dimensions {d = 3). 

We start with a quenched random potential with a Gaussian distribution with variance 
g at each lattice site. The lattice constant, being equal to the monomer size, is a and we 
will measure all lengths in units of a. We will also measure all energies in units of kT which 
is equivalent to setting kT = 1. We first recall the beautiful argument given by Gates and 
Ball 1^ for an ideal chain with no self-avoiding interaction. If we coarse-grain the system 
and divide it into regions of volume ~ R^, the (average) potential in such a region has a 
Gaussian distribution with variance g/R^- Thus the distribution of the potential strength 
V for such coarse-grained volumes is given by 



Here and in what follows we omit constants of order unity. If the volume of the system 
is V R^), the lowest (negative) value U of the coarse-grained potential that is likely to 
be found in such a volume is given by 




(6) 




(7) 



which yields the estimate 



R'uygr^\og{V), 



(8) 



or 




(9) 
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This is the energy per each monomer which resides in this low potential region. Thus if our 
ideal polymer chain is assumed to have a spherical shape of diameter ~ R, it will curl itself 
in the region of the lowest potential given by Eq.(^) above. To determine the volume R^ 
occupied by the chain we recall that we have to pay in confinement entropy for compacting 
the chain below its natural size. Thus we have to balance the confining entropy versus the 
gain in potential energy from the random potential. Thus the free energy is given by 



Notice that an unimportant constant, resulting from the entropy of a free chain, has been 
dropped (5*0 = L log{z) on a lattice with z the number of nearest neighbors, which contributes 
—Sq to the free energy. This term does not depend on R and does not contribute to the 
binding energy). The optimal size Rm is found by minimizing F and is found to be 

~ jki(v) " m- 

where we defined the volume dependent disorder strength by G{V) = g\n(y). Substituting 
this result in F we obtain 

L 



We see that —Fm/L is the binding energy per monomer, and it is strictly positive, so the 
polymer is localized. In what follows we will assume that g is small enough so that G 
for the given system volume, hence Rm ^ 1, and the chain is not totally collapsed unless 
V ^ oo. 

We now add a self avoiding interaction and assume first that it is small, i.e. v « oi 
at least v < g. If the chain is still localized in the same well, which we will see momentarily 
not to hold when L is large, then 




L /^log(V) , VL^ _ r^f.,^2r , „,r^/,n3r2 



.3 +^^-G{VyL + vG{VyL\ (13) 

Here, beside assuming that v is small we assume for the moment that L is not too big so 
the last term in the free energy, resulting from the self-avoiding interaction, is small enough 
so one does not have to take into account the change in Rm due to the presence of f . If we 
plot F vs. L, we see that it is lowest when 



7 



Thus if <^ 1, we have Lm S> 1, and 

For L = 2Lm the free energy vanishes and for larger L it eventually increases fast like 
L^. We can now verify that if L does not exceed Lm then the approximation used above, 
assuming that Rm does not change appreciably from its ideal chain value, is justified. If we 
differentiate the above F in Eq. (|T3|) with respect to Rm we find 



vL 
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again omitting constants of order unity. Thus the correction vGL evaluated at L = Lm is 
of order unity, and we can still use the value Rm ~ l/G. 

For larger L the approximation seems to break down, but fortunately what happens is 
that since the free energy increases when L exceeds Lm-, it is energetically favorable for part 
of the chain to jump into a distant well. Even though there is a cost for the polymer segment 
between the wells one still gains in the overall free energy from the binding energy in the 
wells. Thus the picture that emerges is that as L increases, the chain divides itself into 
separate blobs with connecting segments. In each blob the number of monomer does not 
exceed L^, which is the optimal value for that well. The idea is depicted in Figure 1. 

To be more specific we will now construct a model for the free energy of the chain. 
The first blob will be located in the deepest minimum in the total volume, whose depth 
is roughly given by —G(yY per monomer, where V is the total volume of the system. 
Subsequent blobs will reside in the most favorable well within a range Y , which has to be 
taken self-consistently as the length of the jump. Thus within a range Y the chain is likely 
to find a potential minimum of depth —G{Y^Y. The farther you jump, the deeper well you 
are likely to find. Thus assuming for simplicity that the jumps are roughly of equal size, and 
there are K blobs in addition to the initial blob, the free energy of the chain will be given 
roughly by 

m, Y- L) ^ + + + - -^(^')' + --^'GiY^f) ^ (17) 

with 



K = Lo = G{Y') = 3gHY). (18) 

w + m 2vG[V) 
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FIG. 1: Illustration of how the low free energy conformations evolve with increasing chain length. 
The chain length increases with the direction of the arrows. 

We defined w to be the number of monomers in each blob, and m to be the number of 
monomers in each connecting segment. The term F^/m results from the "stretching" entropy 
of the segment and m/Y^ from the confinement entropy. The term vw? /Y^ represents the 
self- avoiding interaction for the connecting segments. Lq is the number of monomers in the 
initial blob whose free energy was taken care of separately. It is evident that when L is very 
large we can neglect the free energy of the first blob and also take K k, L/{w + m). Thus 
we find for the free energy per monomer 

r/ F{w,m,Y;L) 1 fY'^ m vni^ ^/t^s\2 2^/-,^s\s\ 

L w -\- m \ m Y-^ Y"^ J 

(19) 

This function has to be minimized with respect to w, m and Y to find the parameters 
giving rise to its lowest value. For the connecting pieces of the chain we did not include 
a contribution from the random potential since it is expected to average out to zero for 
these parts. To gain some feeling into the behavior of this function and the values of the 
parameters which minimize it, we display in Table I the value of the parameters and free 
energy per monomer for g — 0.05 and various values of v, as obtained from a minimization 
procedure. The delocalization transition is the point where / changes sign from negative 
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V 


Y 


m 


w 


f 


0.00001 


2206 


2413 


16178 


-0.835249 


0.0001 


346 


534 


2580 


-0.421023 


0.001 


60 


148 


461 


-0.164673 


0.01 


13.7 


66 


97 


-0.0370883 


0.02 


11.1 


69.2 


60.8 


-0.0156444 


0.03 


12.8 


105.5 


41.9 


-0.0056033 


0.04 


22 


302.6 


26.8 


-0.0009992 


0.045 


35.6 


720.1 


20.7 


-0.0001776 


0.0478 


69.4 


2342.4 


16.5 





0.048 


77 


2800 


16 


0.00000725 



TABLE I: Parameters and free energy for the case g=0.05 

to positive, as discussed earlier. Actually, to be more precise, the delocalization transition 
occurs when / = — (lnV)/I/ for finite L, when the translational entropy starts to exceed 
the binding energy. In the hmit of large L we can say that the transition is at / = 0. We 
observe that the delocahzation transition occurs at v = 0.0478 which is close to the value 
of g. We also observe that m <^ w ioi v <^ g and m ^ w near the transition. Also for 
small f , m ~ y, whereas near the transition m ~ Y^. If we compare the value of w from 
Table I with the value of — l/{2vG{Y^)) we find that w is smaller than in the entire 
range. The ratio w/Lm varies from ~ 0.4 to 1 as changes from 10~^ to 0.048. Thus the 
assumption we have made previously concerning this ratio is justified a posteriori. 

Luckily it was possible to solve the minimization equations analytically almost entirely 
in both the limits v <^ g, and near the transition when / ~ 0. Details of the solutions are 
given in the Appendix. Here we only display the results: 

A. The case v <ti g. 
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The parameters are given by 

y = -(lnF-l)-i/2(lnF + 3)-3/2, (20) 

V 

/ = -9/(lnF)2(lnF - l)(lnF + 3)"^ (23) 

The first equation can be easily solved numerically for Y for a given value of v and the result 

substituted in the other equations. Very good agreement is achieved with Table I for small 
values of v. 

B. Solution near the delocalization transition. 
Let us define the parameter 

Av 

(24) 



In terms of this parameter we have 




8^;2 1ny 2vG{Y^)' 
The transition point is obtained by solving the equation 



(25) 

(26) 
(27) 



^ + ^^+ - +1 l-21nP^^+X^^ I I ^0. (28) 

Once the solution Kc is determined for a given g, then the transition point Vc is determined 
from Vc = 3gKc/4:. For g in the range 0.01-0.2 we find that Kc is a number of order unity 
(varies from 1.7 to 0.96 as g changes in that range). This means that is quite close to 
g. Once Vc is known, all the parameters Y,m and w at the transition are determined by 
the solution above. For g — 0.05 we get Kc — 1.21 AA, and thus Vc — 0.0478 in excellent 
agreement with the minimization results from Table I. 

For V > Vc the chain is delocalized. The above expression for the free energy may no 
longer be accurate, but the general picture is clear. There will be very few monomers in the 
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FIG. 2: A typical chain conformation when v ^ g. The dark regions are regions of low average 
potential. Only short segments of the chain are situated in these regions. 



low regions of the potential, and the chain will behave very much like an ordinary chain with 
a self-avoiding interaction in the absence of a random potential. Any little perturbation can 

cause the chain to move to a different location in the medium (see Fig. 2). 

For small values of v when the chain is localized, we still expect its size to grow like that 
of a self-avoiding walk. Thus we expect roughly 

Rg ~ (29) 

since Y is the step size, and K is the number of steps. But since K ^ L/{m-\-w) for large 
L, we find 

Rg ~ 7 ^-T77^^°'''- (30) 

Thus the chain behaves as a self avoiding walk with an effective step size per monomer given 

by 

^ - (31) 

From Table I it becomes clear that the effective monomer size changes from a value of 6 for 
V — 0.00001 to a value of ~ 0.66 at the transition. The reason for the large value of the 
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effective monomer size at very small v is that the chains makes long jumps to take advantage 
of deep wells of the random potential, very much like anchored chains in a random potential 
that make sub-ballistic jumps 0. This is also the reason why the number of monomers w in 
each well is somewhat less than , since a sufficient number of monomers need to be used 
for the connecting segments. For f ^ f c in the delocalized phase we expect the effective 
monomer size to be about 1, since the chain behaves almost like an ordinary self- avoiding 
chain with the random potential not playing any significant role. Thus the chain is expected 
to have its smallest size in the vicinity of the transition. 

It is important to notice that the discussion above is in the limit for very large L. If 
L < l/(2vG{V)), then in the localized phase when f -C f c the polymer will be confined to 
a single well and will appear compact, even though it will not remain so for large L. This 
may explain why in simulations that were done typically with L < 320 the delocalization 
transition appeared as a transition from a compact to a non-compact state of the polymer. 

We should also note that for an annealed random potential there is a transition from a 
collapsed state into an ordinary self-avoiding chain as v increases through the point v = g 
0. This is because the annealed free energy reads 

For V < g the fourth term is negligible and the free energy is lowest when R (when L 
is large). For v > g the first term is negligible, and the radius of gyration grows like 

i?^ (^_(^) 1/32,3/5. (33) 



III. THE CASE OF HARD OBSTACLES 



The case of hard obstacles was investigated recently by us in Ref. [|T^. As already 
discussed in the introduction, three different behaviors were identified as a function of the 
system's volume. Region I is defined when the system's volume V < Vi ~ exp(ci/ \^), where 
Ci is a constant of order unity. Here < x < 1 is the average concentration of obstacles per 
site (total number of obstacles divided by total number of sites). We also assume that x 
is less than the percolation threshold {xc = 0.3116 for a cubic lattice), so sites occupied by 



obstacles don't percolate. In Region I we recall |]14| that in the absence of a self-avoiding 
interaction, the free energy per site for a chain situated in a spherical region of volume R^ 
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in three dimensions is given by 

Fi/L = - \n{z) + + X (34) 

where x is the actual concentration of obstacles in that region, whose minimal expected 
value in a system of total volume V is 



Xm^ X - y ^3 • (35) 

The binding energy per monomer inside the blob resulting from the lower concentration of 
obstacles in this region is given by x — x, since it is equal the entropy gain from a lower con- 
centration as compared to the average (background) concentration x. The chain is "sucked" 
towards regions with low concentration of obstacles since it can maximize its entropy there, 
and these regions of space act like the negative potential regions of the Gaussian random 
potential. Thus in order for the free energy per monomer to reflect correctly the binding 
energy of the chain inside the blob, both the constant x of the background and the constant 
term — ln{z), which is always there regardless of the chain's position, have to be subtracted. 
The relevant free energy per monomer situated in the blob (which is equal to minus the 
binding energy) is given by 



This result coincides with Eq. ([T0|) upon the substitution g ^ x. Thus all the results of the 
previous section carry on to Region I with this simple substitution . 

Therefore we are going to discuss the situation when the system's volume is greater than 
Vi (Region II). In this case the many blob picture still holds, where the blobs are now 
situated in regions free of obstacles (with x = ) whose size is determined again by the 
distance Y of the jump which is also assumed to satisfy ^ Vi ( an assumption which will 
be justified a posteriori). In this case the farther the jump, it is more likely for the chain 
to find a larger space empty of obstacles, which will reduce further its confinement entropy 
w / B? and also the self-avoiding energy ( but there is a cost resulting from the connecting 
segments and the constraint of the total length being fixed). The blob size is given by 



Rmii = ^/Go (the largest expected empty region in a volume ) with JM 
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V 


Y 


m 


w 


f 


0.00001 


552 


2206 


72092 


-0.062 


0.0001 


127 


565 


9135 


-0.051 


0.001 


38 


206 


1241 


-0.034 


0.01 


23 


230 


207 


-0.0077 


0.02 


32 


536 


137 


-0.0019 


0.03 


59 


1737 


120 


-0.00035 


0.04 


194 


13915 


131 


-0.0000099 


0.041 


248 


21327 


134 


-0.0000023 


0.042 


357 


39310 


144 


0.0000024 



TABLE II: Free energy and parameters for a polymer in a sea of blockers when V > Vi. 

Thus the free energy per monomer of a chain consisting of a number of blobs with connecting 
parts, each of length Y, is given by 

w + m\m Y'^ Y-^ J 
where again the constant which is independent of m, w and Y , has been subtracted. This 
assures that the delocalization transition again occurs at / = and not at f — x. The 
results of a numerical minimization of this free energy is displayed in Table II for the case 
of a; = 0.1. 

The transition occurs between v = 0.041 and 0.042. Again we could find analytically 
almost the entire solution both for f <^ a; and at the transition. The solution is given in 
the Appendix. There we show that the transition occurs at v — 0.04142. We observe that 
as for the case of a Gaussian random potential the ratio w/m changes from ^ 1 to ^ 1 
as V approaches the transition from below. The values of Y are seen to be consistent with 
the assumption F'^ > Vi . We also checked that the free energy from Table II is lower than 
what one would obtain by constraining Y to be in Region I, i.e. Y^ <Vi. 

Finally for V > V2, where V2 — exp (^x'^/('^+'^)L^/i'^+^)^ (Region 111), we expect the behavior 
of the chain to stay the same. This is because for jumps within a volume V2 the situation 
reverts to the previously discussed scenario, since the effective volume of interest which 
determines the statistics of the free spaces is of order Y^ and we don't expect Y to be that 
large. 
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An important point to note is that if one performs a simulation with strict self-avoiding- 
walks on a diluted lattice (with x < 1), one has f ~ 1, and hence one will always be in the 
delocalized phase and will not see any localization effects [p!9|] . 

A few words are in order about the case of an annealed potential. This case has been 
already investigated in the literature and we will review it briefly. The free energy in 
the annealed case reads (for d = 3) 

FiR)^^ + xR' + ^^. (39) 

The second term represents the entropy cost of a fluctuation in the density of obstacles that 
creates an appropriate spherical region of diameter ~ R. It was assumed that the chain 
occupies a spherical volume, or at least deviations from a spherical shape are not large 0. 
In the case v = one obtains by minimizing the free energy that R ~ {L/xY^^, a well known 
result. For f > the first term is irrelevant (and so is the "stretching term" of the form 
R^ / L) and one finds that R ~ [v / xY^^L^^'^ . In d dimension, the size scales like L^^'^ when 
f > 0, which is larger than the L^l^'^^'^^ dependence in the f = case. There is no indication 
for a phase transition in these arguments, although some authors speculate that it breaks 
down for large v and a transition to a Flory L^/^''-^'^^ dependence takes place. 



IV. CONCLUSIONS 



In this paper we have considered the case of a long molecule subject to a combination of 
a random environment and a self-avoiding interaction. We have found that in the limit of 
a very long chain and when the self-avoiding interaction is weak, the conformation of the 
chain consists of many blobs with connecting segments. The blobs are situated in regions 
of low potential for a Gaussian random potential or in regions of low density of obstacles 
in the random obstacles case. We have found that as the strength of the self-avoiding 
interaction is increased the chain undergoes a localization-delocalization transition since the 
binding energy per monomer is no longer positive. The chain is then no longer attached to 
a particular location in the medium but can easily wander around due to any perturbation. 
For a localized chain we have estimated quantitatively the expected number of monomers 
in the blobs and in the connecting segments. 
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APPENDIX 

In this Appendix we solve the minimization conditions for the free energy, both for f <^ 
and at or just bellow the transition. The equations obtained by minimizing / from Eq. 



(p!9|) with respect to Y, w and m are in that order 



2Y 2m ?,vm? 6gwG 9vgw'^G'^ . . 

— - y? - -Y~^ — + Y " ' ^ ^ 

-G^ + 2wvG^ = /, (41) 

Y^ 1 2vm ^ 

-— + + ^ = (42) 



m 



where we used 



^^^(F^) = 3^1nr, (43) 

and 

/ = —— — + — + -wG^ + vw^G^ . 44 
w + m \ m Y"^ J 

First, in the case v <^ g, some terms are very small and can be neglected (compare with 

Table I). The equations become 

2Y QgwG 9vgw^G^ ^ ,,,, 

^ + = 0, 45 

m Y Y 

-G^ + 2wvG^ = /, (46) 

-—2 = /' (47) 

f = (^-ujG^ + vw^gA . (48) 

w + m \ m J 



These equations can now be solved analytically with the result given in Eqs. (pO| , pT| , |22| , 

At or very near the transition we put / = on the rhs of equations (^,^). Eq. (^Tj) 
yields 

1 



w 



2vG' 
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The other two equations become 

2m / 1 2vm\ vm^ 3g . 

Thus using the second of these equations to ehminate the first term in the first, the first 
equation becomes 

where we defined k = Av/ {3g) and = v'^m. The second equation becomes 

K^l^x"^ -2x- K^'^ = 0, (52) 

which has the solution 



K 



The values of m and y follows from the relations m = x^/v'^ and Y = k^^'^x'^/v, which 



coincide with equations (p5| , p6| , p7D . Eq. (pH]) just follows from Eq.(Q) and the condition 

/ = o. 

For the case of random obstacles and V > exp(l/ ^/x) we could also solve the minimization 
equations of the free energy both for f ^ x and at the transition. In the first case the 
equation reads 

2Y _ 2wGl _ vw^ ^ 

m syiny Flny ' ^ ' 

-X + GI + 2wvGl = /, (55) 

-— = /, (56) 

f = (— -XW + wGl + vw^GoA . (57) 

w + m \m J 

where 
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We display the solution in the form vahd when v is small enough so that InF ^ 1. It reads 



_1 1 
^"^9(lny)2y^3Gf' ^'""^ 

1 2^, 1 
"^=^9(lny)^.-Gr 

/ = -X + G^. (62) 

Since Go still contains InF, the first equation has to be solved numerically for Y and then 
the rest of the solution follows. 

At the transition the minimization equations read 

— — - + 2wGo—Go + Svw Go— Go = 0, (63) 

m Y"^ Y^ ay ay 

y2 1 2vm ^ 
-^ + y, + ^ = 0' (64) 

-X + GI + 2wvGl = 0, (65) 



with 



Defining 



Axv 



{3x + Gl)ix-Giy 



the solution reads 



va \ a 



2vGl 

The transition point {v^) is determined from the equation 



(67) 



y=^(i±^V. (68) 



- = (70) 



— + 7^ + - ^(^ - ^o) + vw'Gl = 0, (71) 
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which becomes 

l^VTT^ ^ (72) 

Actually since Go and a still depend on InF one has to solve numerically two coupled 
nonlinear equations . The solution gives Vc — 0.04142, y — 302.8, m — 30145, w — 140. 



20 



[1] L. Liu, P. Li, and S.A. Asher, Nature (London) 397, 141 (1999); L. Liu, P. Li, and S.A. Asher, 

J. Am. Chem. See. 121, 4040 (1999). 
[2] D.S. Cannell and F. Rondelez, Macromolecules 13, 1599 (1980). 
[3] G. Guillot, L. Leger and F. Rondelez, Macromolecules 18, 2531 (1985). 
[4] M.T. Bishop, K.H. Langley, and F. Karasz, Phys. Rev. Lett. 57, 1741 (1986). 
[5] A. Baumgartner and M. Muthukumar, J. Chem. Phys. 87, 3082 (1987); see also review chapter 
by these authors in Advances in Chemical Physics (vol. XCIV) Polymeric Systems I. Prigogine 
and S. A. Rice editors, (John Wiley & Sons, Inc., New York, 1996) and references therein. 
[6] S. F. Edwards and M. Muthukumar, J. Chem. Phys. 89, 2435 (1988). 
[7] M. E. Cates and R. C. Ball, J. Phys. (France) 89, 2435 (1988). 
[8] T. Nattermann and W. Renz, Phys. Rev. A 40, 4675 (1989). 
[9] J. D. Huncycutt and D. Thirumalai, J. Chem. Phys. 90, 4542 (1989). 
[10] K. Leung and D. Chandler, J. Chem. Phys. 102 (3), 1405, (1995) ; D. Wu, K. Hui, D. 

Chandler, J. Chem. Phys. 96 (1), 835, (1992). 
[11] J. Dayantis, M.J.M. Abadie, and M.R.L. Abadie Computational and Theoretical Polymer 

Science, Vol. 8, 273 (1998). 
[12] Y. Y. Goldschmidt, Phys. Rev. E 61, 1729 (2000). 
[13] Y. Shiferaw and Y.Y. Goldschmdit, Phys. Rev. E 63, 051803 (2001). 
[14] Y.Y. Goldschmdit and Y. Shiferaw, Eur. Phys. J. B 25, 351 (2002) 

[15] I. M. Lifshits, S.A. Gredeskul, and L.A. Patur, Introduction to the Theory of Disordered 
Systems (Wiley, NY, 1988); I. M. Lifshits, adv. Phys. 13, 483 (1964). 

[16] M. Doi and S. F. Edwards, The Theory of Polymer Dynamics (Oxford University Press, 
Oxford, 1986). 

[17] C. Domb and G. S. Joyce, J. Phys. C 5, 956 (1972). 

[18] P. G. de Gennes, Scaling Concepts in Polymer Physics (Cornell Univ. Press, Ithaca, 1979). 
[19] K. Barat and B. K. Chakrabarti, Phys. Rep. 258, 377 (1995). 



21 



